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Abstract 

We generalize key aspects of gr-qc/10 10.5364 (and also gr-qc/1010.5327) to the case of 
massless quantum field theory on deSitter spacetime. As in that paper, our key objec- 
tive is to derive a suitable "Mellin-Barnes-type" representation of deSitter correlation func- 
tions in a deSitter-invaiiant state, which holds to arbitrary orders in perturbation theory, and 
which incorporates renormalization. The representation is suitable for the study of large dis- 
tance/time properties of correlation functions. It is arrived at via an analytic continuation from 
the corresponding objects on the sphere, and, as in the massive case, relies on the use of graph- 
polynomials and their properties, as well as other tools. However, the perturbation expansion is 
organized somewhat differently in the massless case, due to the well-known subtleties associ- 
ated with the "zero-mode" of the quantum field. In particular, the correlation functions do not 
possess a well-defined limit as the self-coupling constant of the field goes to zero, reflecting 
the well-known non-existence of a deSitter invariant state in the free massless scalar theory. 
We establish that generic correlation functions cannot grow more than polynomially in proper 
time for large time-like separations of the points. Our results thus leave open the possibility of 
quantum induced IR-instabilities of deSitter spacetime on very large time-scales. 
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1 Introduction 

Perturbative calculations of quantum field theory correlators in expanding FRW-spacetimes play 
an important role in cosmology, e.g. for the quantitative understanding of the finer details of the 
formation of density perturbations from quantum fluctuations in the Early Universe. Specifically, 
one is interested in this context mostly in correlation functions of the type 

((^(?,Pl)---(^(?,P£))^ , (1.1) 

where ^ is some quantum field — typically related to the inflaton — and where p,- are spatial mo- 
menta associated with a slice in a flat FRW-spacetime ds^ — —dt^ + a{t)^dx^. ? is a suitably 
chosen time, taken usually when all physical wave numbers Pi/a{t) are considerably smaller in 
magnitude than the Hubble rate d{t)/a{t). For the power spectrum in the CMB, £ = 2 is relevant, 
for "non-Gaussianities" £ = 3 or higher. The state ^ in question is typically fixed as an "in-state" 
in the very Early Universe. It is of considerable interest to: 

1 . Develop methods to calculate such correlation functions systematically in renormalized per- 
turbation theory. 

2. Understand the extent to which the precise choice of the in- state affects the result, and — a 
closely related question — to investigate the IR-behavior of correlators for large times. 

A spacetime which is particularly relevant in inflationary cosmology is deSitter space, with scale 
factor a{t) = e'^ . This spacetime has the additional merit of possessing the same number of 
continuous symmetries as Minkowski spacetime. One may therefore hope to be able to develop 
a formalism for perturbative calculations of comparable efficiency and elegance as in Minkowski 
spacetime. 

In [14J and fi9\, significant progress in these directions was made for the case of massive 
scalar field theories on deSitter spacetime. In these papers "Mellin-Bames"-type representations 
were derived for an arbitrary Feynman graph contributing to an £'-point correlation function in 
a special deSitter invariant state. Such representations express a correlation function in terms of 
complex powers the deSitter invariants formed from the external points in position space. The 
powers are integrated over contours in the complex plane, thus providing an explicit expression 
for the renormalized correlators of the desired type 1). The momentum space correlators can be 
obtained from these by a further Fourier transformation. The representation allows one to see that 
correlation functions have an exponential decay in position space in time-like directions, such as 
e.g. 

((^(Xi)(^(X2))o = 0(e-"'-^^), (1.2) 

where Xi,X2 are deSitter points separated by a large proper time x, and where "0" denotes the 
deSitter invariant state. Similar results were proved in 0141 [T9l also for arbitrary E-point functions, 
to arbitrary orders in perturbation theory, see also [EOl D^- The positive constant is related to to 
the mass of the theory. Estimates of this form imply a quantum version of the "no-hair theorem". 
This states that the correlation functions in essentially any state will approach those of the deSitter 
invariant state at late times. In particular, expectation values of operators in an essentially arbitrary 
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state will approach those of the deSitter invariant state at late times. Hence, this analysis in partic- 
ular answered question 2) above. For previous contributions to this issue, including other claims, 
see e.g. Il22ll0l|2ll26l|29l[30l[3ll[3l. 

Both papers [fT4l[T9l left open the case of a massless theory, such as 



L 



dn, l>0. (1.3) 



The purpose of this paper is to close this gap and to address questions 1), 2) for such theories. 
With respect to 1), we obtain formulae of a similar type for the E-po'mt correlation functions as in 
the massive theory, see eq. (13.431 ). or the related representation given in appendix O The state in 
which these correlation functions are computed is a deSitter invariant state, which is obtained via 
analytic continuation from the sphere. Unlike in the massive theory, this state has no limit as the 
coupling X — 7- 0, as indeed, there is no such state in the free theory [[O- With regard to 2), we find 
that the exponential decay expectedly no longer holds in the massless theory, but we can show that 
the correlation functions grow no faster than a polynomial in Ht. Unfortunately the order of the 
polynomial can depend on the perturbative order. We believe that it should be possible to obtain 
that the correlator actually decays as a polynomial in Hx, and we present in appendix |A] some 
evidence for the two-point function based on the so-called "Kallen-Lehmann representation" in 
deSitter space [4] . Our results imply that, in practice, correlation functions in arbitrary states are 
expected to decay/increase only over cosmolgical time-scales of the order of the Hubble time. 

The feature of massless deSitter scalar fields that complicates our analysis is that-as we have 
mentioned-the underlying free (X = 0) field theory does not have a deSitter invariant state, unlike 
massive theory. In the Euclidean counterpart of the theory on a sphere, this shows up through the 
presence of a zero mode in the massless case. In order to find a perturbation expansion of the 
correlators of a deSitter state in the interacting theory, one can use a modification of the standard 
perturbative setup [25 1 which takes into account the zero mode, as is also often done for perturba- 
tions of massless fields in 2-dimensional fiat space. 

The difference between the massive/mas sless case can already be seen, roughly speaking, from 
the behavior of the 2-point function, W{nP-;Xi,X2), of the deSitter invariant state in the free theory, 
where > is the mass parameter. As we will recall in the paper, this 2-point function is given 
by a hypergeometric function, whose decay for large time-separation is ~ g-cst./fx jj^g positive 
dimensionless constant is of order est. = 0{np-) for — )■ 0, so the decay gets weaker. In the 
free theory, one cannot set = 0. But in an interacting theory (11.31 ). one can set = as long 
as ?i > 0. For an interacting theory, we argue in appendix ^that the 2-point function has the 
"Kallen-Lehmann representation" 

((^(Xi)(^(X2))o = / dM^ p(m2) W{m'^-XuX2) . (1.4) 

[In appendix lAl we also outline an algorithm how to calculate the weight p(M^) in perturbation 
theory.] A massive theory is characterized by the fact that p(M^) has its support contained in 
[m^,oo) for some > 0. By contrast, in the massless theory, the support contains the point 
A/^ = 0. Because the decay of W{M-^;Xi,X2) in time-like directions is only of order e^'^i^ for 



5 



small M^, it is intuitively clear that the full two-point function in the massless theory cannot decay 
exponentially as it does in the massive theory. 

Notations: The dimension of spacetime is D > 2, and is always integer. "EG denotes the set of 
edges of a Feynman graph G, and 'P'G the set of vertices. ^H(z) , '3{z) denote the real and imaginary 
part of a complex number, d/u is the invariant integration measure on the D-sphere S^, resp. real 
deSitter space dSo, depending on the context. Capital letters such as X denote points in resp. 
dSr) (depending on the context), viewed as embedded hypersurfaces in X -X etc. denotes 

the Euclidean resp. Minkowskian inner product in M^+^ depending on the context. For the rest of 
this paper, except in section l34l we will set 

H=l. (1.5) 

However, H can be re-introduced in all of our formulas by simple dimensional analysis: We have 
[X] = 1, [4)] = -D/2 + 1, [m] = -1, [X] = -D + 2n{D/2 - 1), [Z] = 0, [x] = 1, [d/j] = D for the 
length dimensions. Appropriate powers of H then have to be inserted to match the dimensions, 
with [H] = —1. For example, to make the point-pair invariant Z dimensionless [cf. eq. (|3.40l) 1. 
we should set Z = H^Xi ■ X2, the free-field covariance [cf. eq. (IB. 1041) 1 should have dimension 
[C] = 2 [(])], so on the right side we should multiply by H^^, etc. 



2 Massless Euclidean quantum field theory on 
2.1 Euclidean path integral 

Correlation functions in a Euclidean "quantum field theory" on the D-dimensional sphere, S^, are 
formally defined by a path integral 

((^(Xi) . ..<^{Xe))o = ^ /^^ ^(^1) • ■■'^(^e) . (2.6) 

Here / is an action. In our case it is of the form 

1 



^(V(^)2 + P((^) 



djj (2.7) 



with some polynomial potential /'((])) . The path integral has no chance to make sense unless P((|)) > 
est., i.e. unless the interaction is stable. The example to which we will stick from now on is 

75(4)) = + ■ (2.8) 

For > 0, stability requires that ?i > (with 'k = allowed). For < 0, stability requires that 
X > (with X = not allowed). Even if these condition are satisfied, the path integral remains 
only a symbolic expression at best, because the infinite dimensional Lebesgue measure D(|) does 
not exist. It requires a lot of effort to make sense of the path integral, and this has been achieved 
in fact only for D = 2, 3 (and n = 2 in D = 3), and mostly in flat space, see e.g. [9|. In D = 4, the 
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path integral probably does not exist unless X = ("triviality problem"), meaning that the theory 
probably has no UV-completion in the non-perturbative setting. 

At any rate, in practice, one cannot perform the path integral exactly. Instead, it is normally 
interpreted as the formal power series 

In this expression, Iq is the free action, given by 

k = f^^ Q(V(^)' + ^mV) d^i . (2.10) 

While in the exact path integral, the value of np- was arbitrary as long as X > 0, this is now no 
longer so for the individual terms in the perturbation expansion. In this paper, we are interested in 
the massless case, which is 

n? = Q. (2.11) 

This choice creates a problem in the "naive" perturbation expansion (12.91) . as can already be seen 
at the first term V = 0. Then, e.g. for £ = 2, we formally only have to do a Gaussian integration, 
namely (up to a constant) 

j D<^ e-^o ^{X,mX2) = (^-^^ {XuX2) . (2.12) 

Here, is the Laplacian on S^. Unfortunately, this operator has a kemel-the constant functions, 
or "zero modes"-so its inverse actually does not exist! Thus, already the first term in the pertur- 
bation expansion does not exist. This reflects the well-known obstruction to define a "Euclidean 
vacuum state" for a massless scalar field. The same problem also occurs for higher orders in the 
naive perturbation expansion around such a field. 

To see this somewhat more systematically, let us decompose the field <^{X) into the contribution 
from the zero mode called x, and the rest called \\f{X): 

oo N{DX) 

m=x + y\fiX), v|/(X) = £ £ (^LmYUX), (2.13) 

L=l m=l 

where Y^m^L = 0, l,2,...,m = I, ...,N{D,L) are the spherical harmonics on (see appendix |B]) 
satisfying — V^y^„ — L{L + D— 1) Y^m- The Gaussian measure then formally becomes, up to a 
constant 

f = I dx I d\c{^) , (2.14) 

where dVci'V) is the Gaussian measure of "covariance C". By definition 

f . ^ A / ^ fo if n is odd, 

dvc{^\f) U'^^iXj) = \ rfvv^ ; ■ ^^.15) 

J j=i [Upahs {ij)C{Xi,Xj) if WIS even. 
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The covariance C is the propagator for \\f and is given by the "Green's function of — minus the 
zero mode", i.e. 

The problem with the naive perturbation expansion (12.91 ) is now manifest: The integration f^^dx 
implicit in each term in eq. (12.91) is not convergent, because the integrand is a polynomial in x. 
Of course, this problem was not present in the original path integral, where the exponential pro- 
vided a strong damping for large \x\. Thus, to remedy the problem, one has to abolish the naive 
perturbation expansion (12.91 ). However, as suggested by [|25ll . rather than going back to the full 
non-perturbative (not expanded) path integral, which is untractable, one can steer a middle path 
and consider only the ^-integration non-perturbatively, while expanding the rest out. This "semi- 
perturbative" expansion looks in more detail as follows. 
First, we write the potential term as 

Pirn) = + 

Xx^" + Xp{xMX)) ■ (2.17) 
To shorten some formulas, it is then useful to redefine 

X ^ wo\{S'^)-^ X . (2.18) 
We now expand the weighted measure in the path integral as 



V 



= e--'" I 1 (^-±^y(^l,Mj . (2.19) 

We insert this into the path integral, and we also replace ^{Xj) = x + 'y\f{Xj) . The result is 

E V+f 
;=1 ,/■=£■+ 1 

Note that the potential p{x,\f) is polynomial in x (and also V[r). Hence, the i^x-integral can be 
performed explicitly using the standard formula 

r dx e-^" = -X-'l^^'Y^sjln) , (2.21) 
7-cx, n 
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for s odd, and = for even. The Jvc(v) functional integral is then performed (formally) using the 
rule eq. (12.151 ) for Gaussian integration. The covariance ("propagator" of \|/) is C(Xi,X2), rather 
than the non-existent Green's function of V^. We hence get an expansion in terms of Feynman 
graphs. Each vertex in such a graph G corresponds to an interaction term in p{x,\\f), whereas the 
propagators are given by C. 

Carrying out these straightforward manipulations results in the following expansion: 

{(^{Xi)...(^{Xe))o = i E cgX"^Ig{Xi,...,Xe). (2.22) 

^ graphs G 

The sum is over all Feynman graphs G with V intemal vertices, and E extemal vertices labeled 
hy Xi, . . . ,Xe. The valence of the intemal vertices is any number from {2,3, ...,2/7}, whereas the 
valence of the extemal vertices is a number from {0, 1}. The power ac of the coupling constant 
differs from the usual power and is due to the Jjc-integrations (|2.21l) . It no longer has to be 
non-negative, nor integer, but is given by 

-G = ^I^.-|^' (2-23) 

where kj is the incidence number of the j-th vertex of the graph G. cq is a constant depending on 
the combinatorial structure of the graph. It is given by 



CG 



|Aut(G)r^ rr ^^ " (2-24) 



_vol(5^)J ' ' " Y{j{2n-kj 



The second factor comes from the symmetry group of the graph. The third factor comes from the 
Jx-integration (12.211) . and from coefficients in the interaction polynomial p, see eq. (|2.17l) . Iq is a 
Feynman integral. It is given by 

Ig{Xi,...,Xe) = { n [diuiXd] n C(X,(,),X,(,)) , (2.25) 

\i=E+l •'S J (e-EG 

where the product is over all edges £ of the graph G, whose source resp. target vertex is denoted by 
s{t} resp. t{t}. As usual, this expression is still formal because the integral has UV-singularities. 
These must be dealt with using renormalization theory. 

The constant IhC in eq. (|2.22l) is defined so that (l)o = 1. In standard perturbation theory, 
fA^^^ simply removes the "vacuum bubble" diagrams G, and consequently does not have to be 
calculated. However, because the integration over x is non-Gaussian in our situation, lA^^^ now no 
longer just removes the bubble diagrams. It has to be calculated explicitly. In our normalizations. 



2.2 Parametric representation 

We now look in more detail at the Feynman integral Iq, defined formally by (12.251) . To get a 
better understanding of this object, we first give a more useful representation of the covariance C 
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defined previously in (|2.16l ). This representation is derived in appendix El It can be stated most 
straightforwardly if we regard points X G 5^ as unit vectors in ]R^+^ i.e. we think of in the 
following as the embedded hypers urf ace 

S^ = {XeR^+^ : X-X =X^ + ---+Xl = l} . (2.26) 

With this interpretation, we write (Zi -^2)2 = (Xi -X2) ■ {Xi -X2) = 2-2Xi ■X2, etc. The formula 
for C that we will use is: 

r(Y Yi- I 1 r , r(D-i+z)r(-D/2+i-z)r(z) 

+ (-47^)^/2+1/2//^ r(2z+i)r(-z+i/2) [(^1-^2)] . 



Here, the contour c goes around the poles of the integrand at z = —D/2 + 1, — 1,0, and ho = 
~ ^re the harmonic numbers. In the following, we will view c as the union of circles Ck given 

by 

Ck-. t^k + ee^'''', e>0. (2.27) 

We can substitute this formula for each propagator C into our expression (12.251 ) for the Feynman 
integral. To get somewhat cleaner looking expressions in the following, it is convenient to absorb 
the constant term and the pre-factor in C into a redefinition in eq. (12.221 ) of Iq to Iq, of fA^ to fA^, and 
of the combinatorial factor cg in eq. (12:241) to CG given by [28]. Then (nc^)o = ^HT^L^""^ cg h, 
with 

X n r(D-i+zf)r(-D/2+i-z,)r(z,) 

r(2z, + i)r(-z, + i/2) ■ ^ • ^ 

In this expression, z, 7 is a pair of vertices, and for each such pair we have set 

Zij= £ (2.29) 

1. e. Zij is the sum of all the parameters associated with lines that connect a given pair of vertices 
z, J = 1 , . . . , y +£■. If there are no lines in the graph G connecting z, j, then the corresponding factor 
in eq. (|2.28l) is understood to be absent. As usual, we assume that the extemal points X,-, z = 1 , . . . , £ 
are pairwise distinct. 

The integral Ig involves contour integrations, as well as integrals over X, G 5^, z = £ + 1 , . . . , £ + 
V . These 'master integrals' are 

/ V+E r. \ V+£ 

They are absolutely convergent for D = 2 and the configurations Zij G C that we need. But for D > 

2, the integrand is in general too singular at coinciding points X,- — Xj to be integrable. These are 
of course the familiar UV-divergences in perturbative quantum field theory. To treat them, one has 
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to use renormalization theory. This was carried out in [[T4l by a method based on [[T2l [ 
The point is that the integrals over the Xj are well-defined (absolutely convergent) in any dimension 
provided that all ^H(zi/) are sufficiently large. Then, it was shown that the resulting function can 
be analytically continued to other values of the Zij provided that a certain "absence of resonance" 
condition holds between the real parts ^R{zij). The condition states that there should exist no 
integer linear combination E^o ^(zij) £ ^- The subsequent Jz^-integration contours c in Iq 



Ig{Xi,...,Xe) = ldz^iMG{Xi,...,XE;{ze}) 

r(D- 1 +2^)r(-D/2-f- 1 -zi)nzi) 



r(2z, + i)r(-z,+ i/2) 

must then be appropriately chosen so that the absence of resonance condition is satisfied along the 
contour. This can e.g. be achieved by breaking each c up into small circles Ck{t) = fc + e^e^"-" 
around = 0, — 1, —D/2+ 1, with a suitably chosen > for each Jz^-integration. 

2.3 Master integrals 
2.3.1 Example graph 

These issues can be illustrated most easily in the most straightforward example, the "star graph" G 
shown in the following picture. 




In the graph in the picture, the integration variable is Xg, and the extemal legs are Xi, ...,Z7, i.e. 
£■ = 7, y = 1. The corresponding master integral (|2.301 l for this graph G was computed in lfT4l[T9ll 
as: 



'J 



ri-zi)...r{-ZE)r{D+zi + ...+ZE) f<V-/- 

r(Z)/2 + £z,-£w,,) Un-z. + l,w,j)llr{-w,j)( (2.32) 

j^i i j^i i<j V ^ / 

The parameters Zi := Zi[E+i) ^re associated with the lines {i{E + 1)) of the graph G. The contours 
in the expression on the right side run parallel to the imaginary axis in such a manner that the 
left- and right- series of poles of any gamma-function are to the left resp. right of the integration 
contour. Thus, 

• 9^(w, ;) < for any i < j. 
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• ^{-Zi + 'Lj^iWij) > for any /, 

• nZiZ-Lj^,w,j)>D/2. 

These conditions are not mutually compatible for all and indeed the derivation of the 

formula is valid only if these conditions can be satisfied. This is the case e.g. if — £ < < 

for small £ > 0. For general values of Zi E C, we must analytically continue our formula. For 
example, if we want to analytically continue all z, to an open neighborhood of Zi = 0, we must 
move the Zi around some poles of the gamma-functions T{—Zi + Y.j^iWij), and this will give rise 
to corresponding residue. Hence, the resulting formula for Mg{Xi, . . ..Xg) valid for ^(zij) in an 
open neighborhood of will be contain further residue terms in addition to the terms on the right 
side of (|2.32l) such as (e.g. for E — 3): 

Mg(Xi,X2,X3) = ... + (2.33) 

4^«/2+i2.+..+z3 ^ i-x,-x2 y^^^-^^ ^ i-x2-x, y^^^~^^ ^ i-x,.x ,Y^^^^-^^ 

r(-^(zi +Z2 -Z3))r(-^(Z2 +Z3 -Zl))r(-^(Z3 +Z1 -Zl)) 

r(-zi)r(-z2)r(-z3)r(D/2 + i/2+zi+z2+z3) 

These additional residue terms result from moving the Zi around the poles. 



2.3.2 General graph 

For a general graph, one can derive formulae of a similar nature. One possibility is to iterate the 
formula for a single interaction vertex, as pointed out in [|T9l . The resulting formula is derived in 
appendix O Another formula of this type based on the use of graph polynomials and their special 
properties was derived in sec. 5 of [14] 0. To state that formula, it is first necessary to introduce 
some notation. First, we introduce a graph G* whose set of vertices consists of the vertices of the 
original graph G, together with an additional 'virtual' vertex, called *. The edge set of the graph G* 
consists of one edge (z*) connecting vertex i with the vertex *, together with one edge (zj) for each 
factor [{Xi-Xj)^Y'J in eq. (12.281 ). i.e. for each pair of vertices i and j that are connected in G by 
at least one edge. In our formula, there will be integration parameters Wf EC labeled by "forests" 
F within a graph G*. A "forest" is defined to be a subgraph F C G* having the same vertices as 
G, but no loops, and a connected component of a forest is hence a "tree". The forests that we 
consider here have either E trees or E + \ trees. In addition, the former forests have precisely one 
tree connecting an external vertex r E { 1 , . . . , £} with another external vertex s E {1, ...,£} or with 
s — *. An example of a graph G, and a corresponding forest is drawn in the following pictures, 
where E = 4. 

'This paper dealt with a massive field. However, the integrals over Xj in the massless case are of precisely the 
same form as in the massless case-the difference is only in the nature of the subsequent Z£-integrations. 
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Xi 



An associated forest F in connecting X2 with X4 can look like this: 




The integration variables vv = {wf} are not independent, but satisfy the constraint 

Zi 



(2.34) 



for any {ij) such that there is an edge £ G G connection i with j, which we simply write as "(zj) G 
G". We pick a subset of the variables w that are linearly independent from these conditions. 
Our choice is the following. First, we pick, for each {ij) G G, a particular forest F(jj^ having 
the property that Fij fl G = {ij), so that the remaining lines in F(jj^ are all from the set {(/*) : 
i = E -\- l,...,E -\-V}. The corresponding variable wf^...^ is the eliminated via eq. (12.291 ). The 
remaining wf's, i.e. the ones for which F is neither equal to for any {ij), nor equal to 
{ (/*) : z = £■ + !, ...jE + y}, is denoted w. Let us define the meromorphic kernel Kq by 



+ L(o-)eGfej - If3(o-) ) + If h'f) 

r(E(o-)eGzo- - L(/j)^GlF^(o-) 



(2.35) 



nFr(-WF) Y\{ij)eG^{ZF3{ii)WF -Z, 



Yl{iMG^{— +I.Fi>{ij)WF) Y{{ii)^G^{-Z, 



Af 



dWF 

'ini 



All sums or products over F in these expressions by definition exclude the forest F = {{i*) : i = 
E + I, ...,E + V}, and they also exclude the forests whose corresponding integration variables 
have been eliminated via (12.341) . With this notation, our formula, adapted from [14], is: 



Mg{Xu...,Xe) 



Jem 



kg{w,z) n m-Xr-x,)r^ 

\<r<s<E 



(2.36) 
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The complex number a„ (vv) in the exponent is 

a«(w)= £ WF, (2.37) 

F connects r,s 

where we sum over all forests consisting of E disjoint trees, one of which is connecting the vertices 
Xr and Xs, see the above picture for an example. The integrals over w are along a multi-dimensional 
contour C(z) such that y{{wf) = est. for all F and such that the following conditions hold: 

• Dl{wf) < OforallF, 

• ^ [If3(o) ^F-Zij] > for all {ij) e G, 

• E0;)eG^[LF3(.;)>^F-Z,7] < (^+l)/2- 

These conditions ensure that the arguments of the gamma-functions in the numerator of Kq have 
positive real parts, and hence no poles. As with the alternative representation given by eq. (IC.114I ). 
one can show that the vv-integrals are absolutely convergent. This follows essentially because both 
formulas are equivalent up to a change of integration variables. The proof of absolute convergence 
of the alternative form (IC.1141 ) is provided in appendix O 

The first condition on the integration contours is not compatible with the other two if ^{zij) > 
for some (ij) E G. However, we must insert the integral formula (|2.36l ) into eq. (12.281) and perform 
the subsequent integrations over zi along the contour c encircling the poles at = ^^^y^, — 1,0 
to obtain Iqi. This means that, '^{zij) will become non-negative. Thus, we have to analytically 
continue our integral formula (I2.361 l in Zij- What happens is that the w^-contours might have to 
be moved across some poles of the Gamma-functions in the numerator of Kq, and we pick up 
corresponding residue. 

It was shown in lfT4l that our formula for Iq renormalization. More precisely, we showed that 
eq. (|2.30l) can be continued analytically to a meromorphic function of the variables Zij in a subset 
of the complex plane where an "absence of resonance condition" is satisfied, stating that no integer 
linear combination Y^Uij G Z. Eq. (12.361 ) gives an expression for this analytic continuation. 

To obtain Iq, we have to further integrate this expression over G EG along the collection of 
circles q = cq U ... U c_£,/2+i [cf- eq- (12.271 )1. The radii £^ must-and can- be chosen so that the 
absence of resonance condition is satisfied. Thus, our procedure to define Iq in effect involves a 
specific "renormalization scheme", and any other scheme will lead to a different prescription for 
Ig that is connected to the one given by adding finite "counterterms" of the appropriate dimension 
to the action, as described in detail for curved spacetime in [12,1 . The type of possible counterterms 
will as usual depend on whether one has a renormalizable, or non-renormalizable interaction. 

3 Masseless deSitter quantum field theory 
3.1 Analytic continuation 

The deSitter manifold can be defined as the submanifold of {D + 1 ) -dimensional Minkowski space 
IrO+i given by 

dSD = {XeMP^^ : X-X = -Xl + ---+Xl=\} , (3.38) 
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with the induced metric. 




In cosmology, one is mostly interested in the subregion of deSitter spacetime sliced by flat sections 
which is covered by the coordinates G M, x G M^^^ ) defined by 



1 



Xq = sinhf + -eV^ 
2 



^1 



In this region, the metric takes the form 



1 



cosh? eV^ . 

2 



(3.39) 



where dx^ is the Euclidean flat metric on M^^^ The cosmological chart covers the half {X^ +Xq < 
0} of dSo, and its boundary is sometimes called the (a) "cosmological horizon". The cosmological 
horizon is also equal to the boundary 37+ (?^) of the causal future of a point i of J^^. The 
conformal diagram for the cosmological chart is: 



: const. 



: const. 



sections 
Xq = const. 




north pole of 5^ 



In two dimensions, this conformal diagram should be thought of as the hyperboloid which has 
been cut along a vertical line of constant angle. We will use the embedding coordinates X in the 
following, but in our final formulas one can easily go to the coordinates (?,x) in the end. 
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As the sphere S , deSitter spacetime dSo is a constant curvature space, but the metric is 
Lorentzian, rather than Riemannian. Allowing X G C^+^ in the definition (13.381) . one obtains a 
complex manifold, dS'^, called complex deSitter space. The complex deSitter space contains both 
the sphere = dS'^ fl {iM. x MP) (i.e. taking Xq purely imaginary) as well as real deSitter space 
dSn = dS^ n (M X MP) as real submanifolds. In the complex deSitter spacetime, we can introduce 
the C-valued "point-pair-invariant" 

Zi2=Xi-Z2, (3.40) 

where the dot ■ is now the Lorentzian inner product. The analytic continuation of the point-pair 
invariant Z\2 with 0-component of any vector X taken as imaginary, Xq — )■ IXq, is given by the 
Euclidean inner product. Therefore |Zi2| < 1 on the sphere, but on deSitter spacetime \Z\2\ is 
unbounded. The values of Z12 in the real deSitter spacetime are closely related to the causal 
relationship between Xi.Xi E dSo- This is visualized in the following conformal diagram of the 
real deSitter manifold, which indicates the values of Z = Z12 for fixed X2, as Xi varies: 



y+(X2) = future of X2 




Physically, one is interested not in the correlation functions of the field theory on 5 , but on dSr,. 
It is natural to conjecture that the latter can be obtained from the former by analytical continuation 
through the complex deSitter space. But it is certainly not obvious from the outset that such an 
analytic continuation must indeed be possible, nor that it will give a set of correlators on the real 
deSitter manifold with reasonable properties. A set of general conditions on Euclidean correlators 
((|)(Xi) . ..(^{Xe))q which ensure that a reasonable theory on the real deSitter can be obtained by 
analytic continuation was given in [3J, similar in spirit to the "OS -reconstruction theorem" [i23l l9l. 
The key condition is that the Euclidean correlators on satisfy a form of ''reflection positivity". 
Unfortunately, the reconstruction theorem assumes that one has constructed the Euclidean cor- 
relators on non-perturbatively, whereas our construction above was essentially perturbative. 
Therefore, one needs to look at this question more directly by inspecting the analyticity properties 
of Ig, as given by eqs. (|2.28l) together with (|2.36l) . or alternatively by (IC.l 141) . These two equa- 
tions tell us that /^(Xi, . . . ,Xe) is a contour integral over vv,z of an expression whose dependence 
on Xi, . . . ,Xe G enters in the combinations (1 — Z^s)^'', where Z„ is the point-pair invariant. 
Since each (1 —Zrs)^''' is analytic in the cut domain C\ [l,oo), one might expect Iq to be analytic 
in a domain of the form 

=^:={Xi,...,X£eJ5g|Z„. gC\[1,oo) forl<r<5<£}. (3.41) 
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This contains many real deSitter configurations, such as configurations where all points are mutu- 
ally spacelike related to each other. To make this more precise, one has to look at the convergence 
properties of integrals such as (|2.36l ). or (IC.l 141 i. This is discussed briefly at the end of appendixO 
While we do not show analyticity in the entire set S/'e there, we can show analyticity in a large sub- 
domain. In particular, we are able to define Iq in the sense of distributions e.g. in the following 
situations: 

(a) All Xy are mutually spacelike related (Z„. < 1 for all r 7^ 5), or 

(b) For a fixed r, Xr is timelike related to all other points Xs.s^r {Zys > 1 for all r ^ s) and all 
points Xs, s ^ r are pairwise spacelike related (Zst < 1 for all r 7^ 5, t). 

The distributional definition of Ig in an open neighborhood of such points is given by an z£- 
prescription. The correctH /£-prescription is to replace the expression (1 — Xr ■ Xs)^'-' , in for- 
mula (12.361 ) (Euclidean inner product) by ( 1 — ■ + iesrs ) (Lorentzian inner product) . Srs,r<s 
is a sign-function. 



We suspect that Iq is in fact analytic in but this would require a more sophisticated analysis 
than that given in appendixO 

3.2 IR-behavior 

The analytically continued Ig, G any Feynman graph, on real deSitter is thus given by [cf. 



where Z„ are the point-pair invariants in real deSitter space, and where the points Xi,. ..,Xe are 
(for example) configurations in real deSitter space dSo described in (a) or (b) above, and not the 
sphere S^. As we just explained, the «£-prescriptions means that one is dealing with a distribution, 
which in general must be smeared with a suitable test-function first, after which £ > is taken 
to zero. Compared to to the Euclidean expressions eq. (12.361 ) and (|2.22l) . the only difference is 
that I — Xr- Xs (Euclidean inner product) has been replaced by 1 — ■ X^ + iesrs (Lorentzian inner 
product), with Srs the sign function given by (13.421 ). As above, z in eq. (13.431 ) stands for the vector 
consisting of all Z£, i G EG. The z-integral is over the contours z£ E c, introduced above. We obtain 
a completely analogous formula if the alternative representation (|C.l 141) for the master integrals in 

^By "correct", we here mean a prescription that gives rise to correlation functions satisfying the "microlocal spec- 
trum condition" of ||5]. 




(3.42) 



eq. (12361) 1: 




a„(w) 



(3.43) 
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app.Ois used instead. The only difference is that the kernel Kg — Kg and the exponent a,-5 — t- a, ,, 
are modified in the way given in eqs. (IC.l 151 ) resp. (|C.l 16I ). Both representations are expected to 
be equivalent. 

If Xr.Xs are time-like related (so that Z„ > 1), the z£-prescription amounts to putting 

[2( 1 - Z„ + /£5,,)]"" = e""('^-'^)«" (^sinh^ ^) , (3.44) 

with x„. the proper time separating the points. 

In typical applications in cosmology, one is interested in the equal time correlators in the cos- 
mological chart, ((])(?, xi) ■ ■ ■(^{t,XE))o, where all point are hence pairwise spacelike related. To 
obtain an expression for this, one has to simply use the expansion of such a correlation function 
in terms of Feynman integrals, (n<l^)o = 'HT^ Y.'^"^'^ cqIg^ and (13.431) . and substitute the expression 
for the point pair invariant Z„ of (?,x^) and (f,Xi) in the cosmological chart, 

[2(1 - Zrs + /£5„)]"" = e^'*"" |x, - x/"" . (3.45) 

In particular, by taking a further Fourier transform 

/G(f,pi;...;r,P£)= l^n/ j e'P'-i+ - +™/^^^^ (3.46) 

in the spatial variables, eq. (13.431 ) thereby provides an expression for the contribution of the graph 
G to the cosmological observable eq. (|l.ll) . This can be made somewhat more explicit using the 
formula for suitable 5GC\Z/2,Z) even, 

y"rf^-4 Ixp-^e'P" (3.47) 
_ o-2.,^ r(-.-^)r(2. + D-l) , . 

r(-2.) "P" 

Then one sees that, 

f „ ^ s^'HEpO r rr r(D-i+z^)r(-D/2 + i-z,)rfa) 

yGir,Pi,... j,P£j ^^^^ r(2z, + i)r(-z,+ i/2) 

■/ Kg{w,z)Fe( ^,...,%;a) . (3.48) 

Jc{z) V a{t) a{t) J 

Here, a{t) = is the scale factor in the deSitter metric ds^ = —dt^ + a{t)^dx^, a stands for the 
collection {ars{w) , r < s}, and Fe is up to pre-factors an ordinary flat space momentum space 
Feynman integral: 

F£(pi,...,p£; a) (3.49) 
= n --"^2--" r(-a„.-^)r(2a„+D-l) |q„|-2a„-z,.i . 
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The q„-integrations are subject to the usual momentum conservation rule 

Pr = q.r - qrs (3.50) 



s:s<r s:s>r 



for all r = 1 In practice the Feynman momentum space integral Fe has to be defined carefully 

using a suitable analytic continuation prescription in the a^^. We will come back to this issue in 
another work. 

Our main aim in this section is to derive the following result: 

Theorem 1. Let r < E he fixed, and let f{Xi, ...,Xe) be a smooth function of E — I real de- 
Sitter points whose support is compact, and consists of configurations of points (Xi, ...,Xe) 
which are mutually spacelike. Let Xr be a point in real deSitter spacetime which is timelike to each 
point in the support of /. Let us define 

x:=sup{x„ : 5 7^r, (Xi, ...,X£) G supp /} (3.51) 

where x„ denotes the proper time between Xr.Xs. Then we have 

/ /g(Xi,...,X£)/(Xi,...,X„...,X£) <cst.x^ (3.52) 

JXi,...,Xr,...,XE&lSD 

for some A^, and a constant depending only on the graph G, and the function /. 

Remark: 1) Since the correlation functions (n^)o = ^^^Y,^"^cg1g are sums of Feynman 
integrals, we get the same growth estimate also for the correlators. 

2j A look at the proof shows that the expression (|3.52l ) actually has an asymptotic expansion 
~ Lo' cst.„x" for large x. 



Proof: To obtain an estimate, we first have to come to grips with the distributional nature of Iq- As 
we have discussed in the previous subsec. 13. 1[ Ig is defined as the boundary value of an analytic 
function in the neighborhood of the configurations of interest, see item (b) in subsec. 13.11 The 
prototype of distributions of this nature are distributions u{x) on M which are boundary values of a 
holomorphic function u{x-\- iy) that is defined for y > (and e.g. small), and which satisfy a bound 
of the form 

\u{x + iy)\<CQ\y\-'' (3.53) 

for some constant Cq, and for x, say, in a compact set C M. The distributional boundary value 
is defined in more detail as follows. Let zq = xq-\- iya^yo > be fixed, and define, for suitable 
complex integration paths in the upper half plane, 

rz rwM rw2 

v(z) := / dwM I dwM-i--- dwi u{wi) (3.54) 

for z = X + iy,y > 0. Then v satisfies the improved bound \v{x + iy) \ < Ci for a new constant that 
depends linearly on the previous constant Cq. Furthermore, d^v{x + iy) = u{x + iy). Hence, if / is 
a testfunction supported in U, we can define the value of the distribution u{f) by the expression 

«(/):= (-1)^ lim f dxv{x + ie)d^f{x) (3.55) 
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and we have the bound 

|«(/)|<Ci I Jx|a^/W| . (3.56) 

A similar construction is possible for distributions u on U C R", which are boundary values in a 
region U + iV C C'\ where V is some convex cone in M", and the same type of estimates hold. 
Even more generally, an analog of this result holds on complex manifolds if we pass to a local 
chart. 

In the case at hand, the manifold in question is a set of configurations in deSitter spacetime of 
the type described in (b) of subsec. 13.11 By analogy with (|3.53l) . we are looking for a bound of the 
form 

\lG{XuX2 + i£e,...,XE + ieiE-l)e)\<Coe-^ , (3.57) 

where e E is in the future lightcone, and where £ > is small. An estimate of this sort is then, 
by analogy with the 1 -dimensional case, seen to imply a distributional bound of the form 



J Xi,...,Xr,...,XE 

< Ci / , |V^/(Zi,...,Z„...,Z£)|, 

J Xt Xr Xf 



(3.58) 



iXi ,...,Xr,...,XE 

where Ci depends linearly on the previous constant Cq. Thus, our aim is to show that Cq < est. x'^ 
for some A'^, because Ci will then satisfy a similar bound, hence proving the theorem. 

In order to bound Iq in eq. (13.571 ). we substitute the representation (12.311 ) in terms of the master 
integrals, Mq, for which in turn we have the representations (Id 141 i and (IC.l 181 ) [or alternatively 
we could also use eq. (13.431 )1. This gives the estimate (IC.121I ) for the master integral, which can be 
stated as saying that Cq < est. Hs ^^r 1 1 - Z„|^"P^(^") for the analog of the above bound (|3.57l) for 
the master integral. For the point-pair invariant we can use that |Z„| < cst.e^, so we get the bound 
Co < cst.exp(xj^sup9^(v„)). This implies the claim of the theorem if the suprema sup£R;(vr5) 
along the integration paths in (IC.l 181 ) can be chosen negative. Whether this is possible or not 
depends on the values of the complex parameters zg on which the master integrals (IC.l 141 ) depend, 
and we now turn to this question. 

According to our prescription, each z^-integration in (12.311) is broken up into several small 
circles around the points k = — (D — 2) /2, —1,0, cf. eq. (12.271) . Since there is one integration 
contour per integration variable ze, the z = {z^}-integral is a sum 

En/ ' (3.59) 

where n — {nf\. The integrals that we need to look at are therefore 
hiXi , X2 + /ee, X£ + 1 ) e) 

n tt-LG-'\^t-M=^t T{lze + \)T{-zt+l/T) 

■MG{Xi,X2 + ize,...,XE + iz{E-\)e; {ze}) . (3.60) 
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The master integrals Mg are given in turn by eqs. (IC.l 14I ). (|C.l 151 ) where the vv-integral is over a 
multi-dimensional contour w G C(z). This contour is defined so that all variables w in (IC.l 141) run 
parallel to the imaginary axis, and such that any argument of a gamma function in the numerator 
of (IC.l 151) has positive real part. Actually, as we have already discussed, such a contour may be 
defined only as long as all < 0, which is the case if all nf < 0. Then, the contour C{z) 

can be chosen so that when w E C{z), then 9^(vv) < and hence 9^(a„(w)) < 0. In terms of the 
form (IC.l 181 ) of the master integral, this means that we can assume 9^(v„) < along the integration 
contours. Hence we get the desired decay. 

The situation is more complicated when some n^'s are = 0, so that some ze'^ are on the contour 
\ze\ = £e around 0. Then we can get > 0, and consequently it may happen that yi{zij) > 

somewhere. In that case, the contour C(z) is not defined, because one cannot achieve that all 
gamma functions in the numerator of (IC.l 141 ) have arguments with positive real part as well as 
y{{w) < at the same time. Instead, the vv-integral of eq. (|C.l 141 ) is now defined by analytic 
continuation in z. Concretely, this is done by moving some of the vv-contours slightly to the right 
across the poles at of the gamma functions in ^g(vv,z), at the price of a corresponding residue. 
Any integral that we obtain in this way is schematically of the following type: 

where bj.dj E Z, aj, cj, ej E Z". The variable C, E C" is a shorthand for the collection of variables 
z,vv, andthex/s stand for the expressions [(1 — Z„)/2] ineq. (IC.l 141) . The contour ^ is acartesian 
product of (a) small arcs dB)j where are some discs in C, or (b) straight lines dMj parallel to the 
imaginary axis where Mj are some left half -planes in C. The residue of the integrand arise from 
the poles of the gamma-functions and occur if one or more of the linear forms {aj, Q + bj E — Nq 
within Di x ...Hi x .... In terms of the original variables vv,z, this is the case by construction at 
most for vv, z in eq. (|C.l 151) are such that a gamma function in the numerator has a pole. 

To see what type of residues we can get from the integral J, we note a residue formula-type 
integral identity {k > n) for holomorphic /: 

/ ,/ r\Ji^?1} r\^h^ = ^2^')" LkA{P)-d^f{P) , (3.62) 

JdBix...xM„ {{ai,Q +bi)...{{ak,Q +bk) y 

where the sum is over the (discrete) set of P G D i x . . . x D„ such that (aj^.P) +bj^ = ... = (aj^^.P) + 
bj^ = for n linearly independent linear forms. A G Nq is a multi-index of dimension n, which 
is summed subject to the condition that |A| := Ai + ... +A„ < k — n. Furthermore, kA{P) is the 
'winding number' 

1 f M-P)^^...{enX-P)^" ^3_^3^ 

^ Ai!...A„! yaDix...x3D„ {{aiX) +b])...{{akX) +bk) 

and ej the j-th basis vector in C", = d^^^ ...3^^". |A| is the order of the pole. Our integral 
formula can be used (formally) to evaluate 7,noting that the contour ^ is a Cartesian product of 
boundaries of left half-spaces (which can be thought of as infinite disks), and boundaries of discs, 
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and noting that 1 /T{z) is holomorphic. To make the argument rigorous, we should approximate the 
half-spaces by finite discs, and in the process, we will evaluate J as an infinite sum of residues at 
points P. The convergence of this infinite sum does not have to be considered in practice, because 
only finitely many residues give the dominant contribution in the large-distance analysis of (13. 601 ) . 

In our case, P consists of all possible vectors with entries (z, w) , such that a linearly independent 
set of the linear conditions is satisfied which state that (z, w) is at pole of a set of gamma functions 
in the numerator of eq. (IC.l 151 ). The number of those conditions has to be the same as the number 
of entries. By construction, any such P will have components with real part < 0. Hence the term 
with the strongest growth as |Z„| — oo in eq. (IC.l 141) is one corresponding to a residue at P with 
all a„-,Z£ = 0. This will produce a term of the form cst.(log(Z„ — 1))^, where is the order of the 
pole. Because log(Z„ — 1) < cst.x, his demonstrates the claim of the theorem. □ 



3.3 Resummation 

An unsatisfactory aspect of thm. \T\is that, according to the theorem, the ii-point functions could 
still grow polynomially in x for large time-like separation x. It is conceivable that better bounds 
could be obtained using more refined methods. For the 2-point function, the methods based on 
spectral representation outlined in appendix |A]might be one possibility. The best option would of 
course be to give a full non-perturbative analysis of the correlation functions, but this seems to be 
very difficult. A more modest option could be to perform resummations of certain infinite classes 
of Feynman diagrams. Let us outline this here at a simple example. The action is rewritten as 



^(V\|/)^ - ^m^(x)\|/^ - (higher order in \|/) 



dfi , (3.64) 



with the usual decomposition of (|) = jc + \|/ into the zero mode and the rest. The quadratic term 
has a "mass" given by nf- = 2n{2n — l)'kjc^ > 0. The simplest possible resummation is to perform 
in closed form the perturbation series generated by this quadratic term. As usual, carrying out the 
corresponding geometric series is equivalent to absorbing the mass-term into the covariance of the 
Gaussian measure. The new covariance is 

oo 

C'W =Co £(-m2)^(Co...oC) (3.65) 
v=o 

From the definition of C in terms of spherical harmonics it then follows 

C'iXuX2;x) = ( \ ^ ) (Xi,X2) - — • (3.66) 

Note that this covariance depends upon x via = 2n(2n— l)'kx^ > 0, and also note that the 
difference on the right hand side gives C'{x) — )■ C as x — )■ 0, and is therefore well-defined. Hence, 
the path integral becomes 

I e-' n^(^./) = f^y^ e"^'" / exp (^-l l^^p>,{x,y\f)d}?j fl[x + y\f{Xj)] , 

(3.67) 
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where p>3 represent the terms in p that are higher than quadratic in \|/. The last exponential 
is expanded in a Taylor series, as usual. This expansion differs from the original perturbation 
expansion, because each term already involves an infinite sum of certain diagrams. 

For E = 2 points, this partially resummed path integral seems to behave better for large time- 
like separation, at least to low expansion orders of the exponential. For example, the contribution to 
the above path integral from the lowest order term in the expansion of the exponential exp(. . . ) = 
1 + . . . is seen to behave as est. loglogZ for large Z. This behavior is better than that given by 
thm.[Tl and it shows that the large Z-behavior might be improved by a partial resummation. Stated 
differently, the unsatisfactory growing nature of the bounds in thm.[T] could be an artefact due to 
the truncation of the perturbation series at finite order. 

A more ambitious program is to try a resummation of larger classes of diagrams, via the so- 
called "skeleton expansion", as considered in [|2U, and thereby to obtain even better bounds. How- 
ever, when performing such expansions, one has to be careful about the correct renormalization 
prescription, as we now briefly explain. The point is that any change in the renormalization that 
we have adopted here, can be absorbed into the addition / — )■ / -f 5/ of finite counterterms to the 
original action, by the general theorems [|5] [l2] [131. The precise form of the counterterms is dic- 
tated by power counting and covariance. Let us focus on the case D = 2 for simplicity, where even 
a non-perturbative existence proof of the deSitter correlators is available [16]. The counterterms 
in 5/ take the same form as the terms already present in the action. In particular, we may get 
a counterterm of the form 5m^ (|)^, which is of course of the form of a mass term. Such a term 
will in effect make the theory massive if "all diagrams are summed" i.e. in the non-perturbatively 
defined theory. It could lead to an exponential decay of the correlators. Of course, this is a radical 
change of the nature of the theory, and it should be imposed that the theory remains massless. In 
deSitter space, this "renormalization condition" is that the support of the Kallen-Lehmann measure 
p(M^) dM^ (cf. appendix lAl) should contain the point = 0. Such a renormalization condition 
must also be respected by any resummation taking into account only a limited class of diagrams. 

3.4 Physical consequences 

For a massive interacting field, the deSitter correlators {<^{Xi)...^{Xe))o in the deSitter invariant 
("Euclidean") vacuum state decay exponentially in time lfT4l [TSl [T9l to all orders in perturbation 
theory, but for a massless field, we have seen that they grow polynomially in time. Although 
our analysis was only graph-by-graph in a semi-perturbative setup, such a behavior, if true non- 
perturbatively, can potentially have significant physical implications for the evolution of the uni- 
verse on large time-scales. To see this, let O be the operator 



where the sum is finite and where ^„ are some "wave-packets". For technical reasons-since we 
want to apply thm. \\\-we assume that the support of each ^„ is compact, and consists of config- 
urations (Xi, of points which are mutually spacelike to each other. Let ( . )>{/ be the state 




(3.68) 
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obtained by applying O to the deSitter invariant state ( . )o, i-C- 

WX.)...^W).:^ <°-^'^;^^^'°>° , (3.69) 

or in "vector notation" (i.e. in the GNS-representation of the deSitter invariant state), |^) = 
C'|0)/||O|0)||i/2. Let y(x) be a time-like curve parameterized by proper time x, which goes to 
future infinity y^, and which eventually becomes timelike related to any deSitter point X, in the 
support of ...,X„) for all n. Then the growth of the deSitter correlators in ( . )o stated in 

thm.[T]and the following remarks immediately give: 

((^[y(x)])>P ~ P{Hx) = cst.o + cst.i (//x) + cst.2(//x)2 + . . . , (3.70) 

where we have reintroduced the Hubble constant H. P is a polynomial, which in our analysis, 
depends on the order to which the perturbation expansion is carried out. The coefficients of the 
polynomial depend upon the precise choice of the wave packets, the value of H, and X. As a 
function of X, the constant term is of order one. The higher terms in the polynomial are of order at 
least X^/^" for a interaction. 

Although we have not analyzed composite operators in this paper, this can be done. One obtains 
e.g. that, to lowest order in X 

(p[Y(x)])vp~(p[y(x)])o + cst.^(//x)2 + ... , (3.71) 

where p = T^js/'^'f is the energy density operator evaluated along the curve, and where est. is a 
constant depending on H and the precise form of the wave packets. For fixed wave packets, and 
small H, the constant would be of order in D = 4. The expectation value on the right side in 
the deSitter invariant state depends on the renormalization convention for the composite operator 
r^v Given that deSitter spacetime ought to be a solution to the semi-classical Einstein equations, 
it is natural to fix the renormalization convention by requiring that 87rG(71jv)o = — A^^v In D = 4, 
we have A = 3//^, so we get that to lowest order in X 



3 „ fH\ 



2 



(3.72) 



with Ep the Planck energy, and est. a constant of order unity. Thus, a perturbative analysis sug- 
gests that the expected energy density for a self-interacting, massless field in a "typical state" (not 
equal to the deSitter invariant state) will grow in time, and could thereby give rise in principle to 
significant back-reaction effects. However the linearly growing term will be comparable to the 
vacuum energy term only when X{H /Ep)^{Hx)^ is of order one, which is satisfied only when the 
time-scale x is of the order of the Hubble time, and when X > (Ep/H)^. Inserting the presently 
observed value for H, this would correspond to a huge value of the coupling constant ?i ^ 1 for 
which our semi-perturbative analysis is clearly not applicable, and a non-perturbative analysis will 
be required to settle the issue. 
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A Spherical harmonics, Gegenbauer polynomials, Kallen- 
Lehmann representation 

Here we outline how one can obtain a "spectral representation" of the 2-point correlation function 
of an interacting field analogous to the Kallen-Lehmann representation in Minkowski spacetime. 
The derivation of this formula involves spherical harmonics in D-dimensions, so we briefly recall 
their basic properties. For more details, see e.g. [2J. For simplicity, H = \ in this subsection. 



A.l Spherical harmonics and Gegenbauer functions 

Spherical harmonics on the unit can be introduced via harmonic polynomials in the em- 
bedding space M^+^. A polynomial P{X) on M'^+i is called homogeneous of degree h if 
P{XX) = X^P{X), and it is called harmonic if it is a solution to the Laplace equation on M^+^. 
The harmonic polynomials of degree h = L form a vector space, the dimension can be seen to be 
N{L,D) = ^)' _ Spherical harmonics on of order L are by definition just an or- 

thonormal basis of the space of harmonic polynomials, restricted to S^. The spherical harmonics 
Yij{X)J= l,...,N{D,L) dst thus normalized so that 

Y,YLm{XlYYLm{X2) = 6(Xl,X2) , / d^iiX) Yt„r{XyYa,^,{X) = huL>K,m' (A.73) 

where the 5 function is that on 5^, defined with respect to the measure di-i. Expressing the Lapla- 
cian on M^+^ in polar coordinates, on sees that the spherical harmonics are eigenfunctions of the 
Laplacian on the D-sphere with eigenvalues —L{L+D — 1), so that L may be viewed as the ana- 
log of the total angular momentum-, and m may be viewed as the analog of the magnetic quantum 
numbers. One has 

N(D,L) or I n_i 

L YUXiYYUXi) = ~Z,oJ , (A.74) 

m=l ) 

where are the Gegenbauer polynomials, and where Z is the point pair invariant. The Gegenbauer 
polynomials are expressible in terms of a hypergeometric function, 

Eq. (IA.74I) may be viewed as saying that the Gegenbauer polynomials are, up to normalization, 
the integral kernels of the projector Ei onto the eigenspace for the eigenvalue —L(L + D— 1) of 
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the Laplacian on S^. Since the dimension of this eigenspace is equal to N{D,L), one gets the 
orthogonality relation \x{EiEii) = N{D,L)dix'- Writing out the trace of the integral kernels as 
integrals, one infers from this that 

l\dZ (1 -Z2)«/2-i ci''-'^/\z)c[^-'^/\z)=No,L Kv (A.76) 

for L,L' G No, with normalization factor 

vol(5^-i)r(L + D-l) 
vol(5^)(2L + D-l)r(D)r(L+l) ' ^ '^^^ 

A.2 Kallen-Lehmann measure p 

These formulas can be used to first obtain expressions for the Euclidean two-point function for a 
field of mass > 0, as follows. Writing out the condition that ((|)(Xi)(|)(X2))o is the Euclidean 
Green's function of (— + m^) on the sphere gives using (|A.74I) 



- voKS^-MlV"- ^ ^-c(c + D-1)+L(L-D+1) ^^-"^^ 

where c = -(D- l)/2+ [(D- 1)2/4 -m2] 1/2. Using c[^'^^^^{Z) = (-1)^ c[^"^^/^(-Z), and 
using the above representation of the Gegenbauer polynomials as hypergeometric functions, the 
above sum can be converted to a contour integral over L with the help of a Watson-Sommerfeld 
transformation, as observed in [fifSl : 

((^(Xi)(^(X2))o = / — (2L + D- 1) Pt ^L{Z) (A.79) 
Jc 2111 

where the contour C is running parallel to the imaginary axis, leaving the poles in the denominator 
of 

m"^ +L[L + D— 1) 
to the left, and the poles at No of to the right, see the figure. 
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C" 



C 
-> 



m{L) 



-D+l-c 



The kernel in this formula is defined as 



Ar.(Z) 



1 r(L + D-l)r(-L) 
(47r)0/2 r(D/2) 



2F1 -L,L + D-l;D/2: 



1+Z 



(A.81) 



We can deform the contour C to the contour C" by moving it across the pole alH L = c. The 
integrand as well as the contour C" is anti-symmetric under L — )■ —L + D— 1, so we are left with 
the residue: 



((^(Xi)(^(X2))o = Ac{Xi ■X2) {= W{m^;XuX2)) , 



(A.82) 



giving the free Euclidean Green's function for the mass parameter = c{c + D — 1). For- 
mula (IA.79I ) can be verified as follows. The function has poles at L = 0, 1,2, .... Then, de- 
forming the contour C to a contour C' that encircles these poles along the positive real axis (see 
figure), we can evaluate the integral by means of the residue theorem. The residue at L = 0, 1 , 2, . . . 
of {2L + D— \ ) Pi Ai is precisely equal to the L-th term in the sum of eq. (IA.78I ). 

Pl is interpreted as the "power spectrum", or "spectral density". For L E Nq, the power spec- 
trum can be obtained by multiplying both sides of (|A.78I) with c'f ^^^'^{Xi ■X2) and integrating 
over Z = Xi ■X2 using the orthogonality of the Gegenbauer functions (|A.761 ). Using also (IA.81I) 
and (IA.75I) . one obtains the inversion formulas: 



r(D) 



r{-L)r{L+D-i)JsD>csD 

{<^{X,mX2))o 



((^(Xi)(^(X2))oAz,(-Xi-X2) 



^ sinTrL , 

y (2L-fD- 1) Pl Al(-Xi ■X2) . 

L=0 ^ 



(A.83) 
(A. 84) 



These formulas are valid not only for the free field correlators, but in fact also for the interacting 
correlators (denoted in this section by {^{Xi)^{X2))ox to distinguish them from the free ones). 



^Here we assume that < {D— 1)^/4, i.e. that we are in the complementary series. For the principal series, a 
similar contour symmetric under L — L + D — 1 is chosen. 
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because in their derivation only the invariance property under 0{D+ 1) was used. The spectral 
density in the interacting theory is denoted P^ x'^ of course it is no longer given by the same for- 
mula (IA.80I) as in the free theory. Thus, we have the inversion formulas 

= r(-L)r (fio- 1) i.,,<^<^-'^<^^))°-^ ^^<-^' -^^^ <^-«« 

°° sinjrL 

(^(Xi)^(X2))o,x = E — — (2L + Z)- 1) Pi^^x Ad-X, ■X2) . (AM) 

The first equation now has to be interpreted as the definition of P^ x. Note that these formulas 
define P^ i not only for natural numbers L, but even provide an analytic continuation for complex 
values of L. The poles in L of Pi^ x give information about the mass in the interacting theory. For 
X = 0, i.e. in the free theory, there are poles at L = c or L = —D + 1 — c. For a principal series 
scalar field, we have c = — (D — l)/2 + zp, p G M, while for a field in the complementary series we 
have c e [— (D — 1 ) /2, 0] . In an interacting theory, the singularities of P^ x lie in these regions, too, 
but they do not have to be simple poles. 

Let us suppose that we knew that there were no singularities in the region ^H(L) > 0, and let us 
also suppose we knew that \Pix\ grows slower than e'^l-'^^^l for large |L|. Then, because the large 
I J(L)| asymptotics of A^ is given by 0{e^^^^^^^^), we can apply the same argument as in the free 
theory to get the representation 

{<^{Xi)<^{X2))o,x = j^^. {2L + D-l)P^^xMZ) (A.87) 

also in the interacting theory. A major advantage of this formula is that one can directly analytically 
continue both sides to real deSitter space; the only difference is that on the right side, we then have 
to take Z — )■ Z + ies, where s is the sign that indicates whether Xi is to the future of X2, or vice 
versa. 

Eq. (IA.87I) may be used under certain conditions to derive an analog of the Kallen-Lehmann 
representation in Minkowski spacetime, because Ai is equal to the free field two-point function 
with = —L{L + D— 1). Indeed, consider first the case that the P^x is analytic all the way to 
9^(L) > — (Z) — l)/2. Then C can be moved all the way to the left to a contour C' parallel to the 
imaginary axis with ^H(L) = — (Z) — l)/2 + e for any e > 0, see the following figure. 
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^{L) = -^ 

The red zigzag line indicates the possible location of singularities such as poles. Suppose P^ x 
has a suitably distributional limit for £ — 0. Then using ^ = P-l+d-iX^ ^^'^ changing the inte- 
gration variable from L = - (D - 1 ) /2 + £ + / [M^ - (D - 1 ) 2/4] to M^, the above formula (IA.87I ) 
becomes 

f'OO 

{^{Xi)^{X2))q,x = L ,,2 dM^ Pp(m2) Wp{M^;X,,X2) , (A.88) 

where Wp{M^ ,X\,X2) is the two-point function on deSitter spacetime of a free field in the principal 
series for mass >{D— 1)^/4. The "Kallen-Lehmann weight" is given as the discontinuity 

across the zigzag line. In the general case, one expects that Pi x can have further singularities for 
L G [— (D — l)/2,0]. Suppose that there are no singularities in a neighborhood of 0. Then, we may 
move the integration contour C to the contour C' depicted in the figure, where red zigzag line again 
indicates the possible location of singularities. 
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The horizontal piece of C' can be parameterized by L = — (D— l)/2±z£+ [{D — 1)^/4 — 
M^] Then, changing integration variable for this horizontal piece to < < {D— 1)^/4 and 
assuming that the discontinuity across the zigzag line 

exists as a distribution, the Kallen-Lehmann representation in the general case has an additional 
piece and is given by 

(D-l)2 

poo n — 

J — — Jo 

(A.91) 

Wc{M^ ■,Xi,X2) is the two-point function of a free field of mass < (Z) — 1)^/4 in the comple- 
mentary series. The first integral comes from the vertical part of the contour C' , and the second one 
from the horizontal part of the contour. The first integral is the contribution to the Kallen-Lehmann 
weight from the principal series, the second is that from the complementary series. Formulas of 
this kind have previously also been derived in the pioneering work [4J. For a free massive field of 
mass > and Pl given by eq. (|A.80I) . we get, using the standard discontinuity formula 

— = Inihix) (A.92) 

x — iz x + ie 

that pp,c (M-^) — 5(M^ — m^) , depending on whether is in the principal- or complementary series, 
as is of course required. 



A.3 General consequences of Kallen-Lehmann representation 

If Pi^ x has no singularities in a neighborhood of 0, then by definition, the Kallen-Lehmann weight 
pc will have its support in supppc C [m^, (D — 1)^/2] for some > 0. It makes sense to think 
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of such theories as "massive", with the value of the mass parameter. On the other hand, if the 
singularities of P^x go the way to L = 0, one would speak of a "massless theory". In that case, 
the above contour distortion argument does not work as stated, as one cannot thread a contour 
C' between and the singular red zigzag line. However, one still expects the formula (IA.911 ) to 
hold true, because it simply expresses the fact that the unitary representations of SO(D, 1) can 
be decomposed into a direct integral from representations of the principal-, complementary- and 
discrete series [27]. Indeed, if the theory has been constructed non-perturbatively, then the positive 
definite property of the 2-point function, 

/ {<\>{Xi)^iX2))o,xf{Xi)fiX2j>0 (A.93) 

can be thought of as defining a positive scalar product on the space C^((iS£)) of compact support, 
on which the group SO{D, 1) acts in a natural fashion. The representation in question is obtained 
from a suitable completion. The discrete series is absent because it is in conflict with the locality 
property required for the 2-point function. The connection to the representation theory also estab- 
lishes that Pc,p(M^) are not only distributions of positive type, but that Pc,p{M^) dM^ are actually 
(positive) measures [14J. 

Let us write p = \{M < ^^j^} Pc + ^{M > Pp for the combined measure on M+ which 
we assume exists, and a be the supremum of all numbers such that 

limsup (r''-^ f dM^ piM^)] < °° ■ (A.94) 

a gives a measure of the decay of the Kallen-Lehmann weight at = 0. For a massive theory, 
the Kallen-Lehmann weight is supported away from 0, suppp C [m^,oo), so a = oo. For a < 0, the 
representation eq. (IA.911 ) would be ill-defined as fW{M^;Xi,X2)f{Xi)f{X2) ~ 0{M^^) for smaU 
M^. For a massless theory, the topic studied in this paper, a is by definition a finite number > 0. 

Proposition 1. Suppose that a non-perturbatively defined scalar field theory has a two-point func- 
tion with Kallen-Lehmann representation. Suppose that Kallen-Lehmann weight has a > 0. Then, 
for large time-like separation (Z oo), the two-point function decays as 

mXimX2))o,x\ < est. (logZ)-« . (A.95) 

Proof: Of course, we use the representation Jq dM^ PiM^) W{M'^;Xi,X2) for the 2-point 
function, where V7(A/^;Xi ,^2) is the 2-point function in the free theory of mass M^, see eq. (IA.82I ). 
We split up the integral into Jq' = /q + J^. For > £ > 0, estimates on the hypergeometric 
function (see e.g. IfTTl ) give the uniform bound 

\W{M^;Xi,X2)\ < est. (A.96) 

for Z — 7- 00. Consequently, the contribution from the second integral is bounded by est. < 
est. (log Z)^" for large Z. For < e, we have instead the uniform bound 

\W{MhXuX2)\ < ^ Z'^' , \W'{M^-Xi,X2)\ < ^ Z-^' (A.97) 
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for Z — 7- oo, where a prime denotes a derivative w.r.t. M^. Therefore, using a partial integration 
together with the assumption (|A.94I) . 



dM^ P(m2) W{M^-XuX2) 



W{M^;XuX2) d 



d// pi/J^ 



r£ ( rM^ 



< est. / JM2m2(«-i)z-^'. 



(A.98) 



The last integral is given by est.(logZ)^'^. This eompletes the proof. □ 
The theorem gives no information when a = 0. However, in that ease we ean e.g define a 
refinement of this quantity by taking a to be the supremum over all numbers sueh that 



limsup f (log?)"^""?"^ [* dM^ p{M^)] < 
t^o V Jo J 



(A.99) 



Sueh a eondition, if satisfied for a > 0, would give a deeay | ((|)(Xi)(|)(X2))o.xl < est.(loglogZ)^'^. 
The proof would be similar, and it is also elear that one ean eonsider many further refinements 
along those lines. 



A.4 Kallen Lehmann measure in perturbation theory 

In perturbation theory, Pc.p is only defined in perturbation theory, and we eonsequently eannot 
apply this result direetly. However, if we knew for example that pc was absolutely eontinuous near 
0, then a qualitatively similar definition of a quantity a as above ean be made, sueh as e.g. the 
supremum over all numbers sueh that limsupy^^2^QM^^°'pc(M^) < oo, and one ean derive a similar 
deeay result. 

Also in perturbation theory, we ean ealeulate Pc.p for the interaeting theory as the diseontinuity 
of Pi x along the zigzag line [ef. eqs. (IA.90I ), (lA. 891 )1, where Pi^ ^ is now ealeulable explieitly graph- 
by-graph. Indeed, breaking up the 2-point funetion in the first equation (IA.85I) into individual 
Feynman diagrams, we have, with Iq the Feynman integral (12.311) . and with 

Plg-=—, r/ /g(^i,^2) Al(-Xi ■X2) (A.lOO) 

r(-L)r(L-fD-i) A^xs^ 2; Lv 1 2; V ) 

that 

Pl,x = Y.^g'>^''''Pl,g. (A.101) 

G 

Consequently, 

(A.102) 



32 



thus giving the desired breakup of pc- into contributions from individual graphs G. A similar 
expression holds for pp. 

Note that the integral in eq. (|A.100I) simply produces a new Feynman integral P^ q = Ig{l) 
where the two external points Xi.Xj of G have been closed off with a new type of propagator, 
A/,(— Xi ■X2), with "mass" = —L{L + D— 1). This gives gives a new graph G(L) without exter- 
nal legs, with a distinguished line carrying the new propagator Ai instead of C, see the following 
picture. 




The new propagator can be represented in Mellin-Bames form as 

^'^-'''-'''^ = i-nr^+yd-e-.^ r(2z + D)r(-z-D/2 + l/2) [(^^"^^^ ^ ' 

(A. 103) 

which has a similar structure as C, see thm. |3l Therefore, in conjunction with the Mellin-Bames 
representations for Iq derived in the body of the paper, this will give Pl,g as /^(l) in the form of a 
Mellin-Bames integral similar to eq. (|2.28l ). One can use such representations to analyze in more 
detail the power spectrum Pn hence pc,p, and to thereby analyze decay of the 2-point function. 
We will not do this here. 



B Expressions for the covariance C 

In this section we provide the representations for the covariance C{Xi,X2) on the sphere used 
in the main text. The discussion is somewhat different for even and odd D; we focus on the case 
when D is even, as in the main text. The definition of C{Xi,X2) is 

C(Zi,Z2) = ^ L(L + D-1) 

Here Yi^^^ = 1, ■■■,N{D,L) are the spherical harmonics on the unit S^, and in the second line 
we used the definition of the Gegenbauer functions cj^ ^^^^{Z), cf. eq. (IA.74I) . We also use the 
shorthand Z — Xi-X2 as in the main text. The above sum can be converted to a contour integral 
over L with the help of a Watson-Sommerfeld transformation, as already explained in a similar 
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context in the previous section: 



C(X,.X.) ^ c..r;|.c.g(.i)i^±§^cr')/^(Z) (B..05) 



-;oo+£ ^iii i^yi^ -r 1-^ — 1-) 



= est. 



jp2pr(-p + ^)r(p+ 2 



2 . 



K2ni (_p + ^)(p + 

X2F,(— + P. — -P;^;— j . (B.106) 

Here e > is small. The poles of ctg(7rL) at L G Z with residue +n generate the original series if 
we perform the first integral by the residue theorem, moving the contour to a contour encircling the 
poles of the integrand at L G Z+. In the second line, we have used eq. (|A.75I) . the transformation 
of the Gegenbauer functions under Z — — Z, the standard identity sin(7u) = 7T:r(l +x)r{—x), and 
we have switched to the integration variable p = L+(D— l)/2. The contour K is running parallel 
to the imaginary axis, with < 9^(p) < ^^j^- 

We note the integrand in the last expression is anti-symmetric under p — )■ —p. Therefore, the in- 
tegral would vanish if the integration contour was invariant under this transformation. The contour 
K is actually not invariant under such a transformation, but we can deform it to such a contour K', 
e.g. taking K' to go along the imaginary p-axis. When we deform the contour K to K', we cross 
the pole at p = (D — l)/2 of the integrand, and we pick up a corresponding residue. To calculate 
that residue, we may use the standard power series expansion for 2F1 : 

r(-p + ^)r(p + V) ^fD-l_D-l ^,Dl+Z 



r{D/2) ^^^V 2 2 P'2' 2 

-„?0 niriD/2 + n) (,^J " ^^''''^ 

As is clear from this expression, the n = 0-term has a pole at p = (D — 1 ) /2, giving rise to a double 
pole in the integral eq. (|B.105I) . whereas the other n > terms are analytic there, giving rise to 
single poles in that integral. Thus, using the residue theorem and keeping track of the constants, 
we get: 

where \|/ = T' /T is the psi-function. We can now perform the infinite sum using the series 
L«>0^"/" = — log(l — x), and we may also use the standard identity \\f{n) —\\r{l) = h„ for the 
psi-functions, where 

hn = '^ + l + l + --- + ^ (B.109) 
2 3 n — l 

are the harmonic numbers. Then we easily get the following result, which is our first alternative 
representation of C: 
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Theorem 2. The covariance C defined by eq. (IB. 1041) is given by 

1 \TiD-l)hD 



C{XuX2) 



(47r)^/2 1 r(D/2) 



(B.llO) 



az 



1 -z 



We can carry out the differentiations in the last expression. If this is done, we find 



C(Xi,X2) 



1 



l_Z\-^/2+iD/2-2 Y{D/l + n) f\-Z 



(471)^/2 t V 2 

Tin -I 



+ 



r(D/2) 



„tj, n\{D/2-\-n 

1 -Z 



+ 



/zo + \|/(Z)-l)-v(Z)/2)+log 



(B.lll) 



We claim that the sum on the right side can again be expressed via a contour integral, which is our 
second alternative representation: 

Theorem 3. The covariance C defined by eq. (IB. 1041) is given by the following contour integral: 



C{XuX2) 



h 



D 



+ ■ 



i fdz 

)Z)/2+i/2 y/^ 



T{D-\ +z)r(-D/2+ 1 -z)T{z) 



[{X1-X2) 



2iz 



voi(5^) (-471)^/2+1/2 y,- r(2z+i)r(-z+i/2) 

(B.112) 

The contour c is encircling the poles at z = —D/2 + 1, —D/2 + 2, ...,0, and ho are the harmonic 
numbers. 

Proof: Using the definition Z = Xi ■X2, the doubling identity T{x)T{x+ 1 /2) = 2-2i^+i ^T{2x) 
as well as r(l/2 + ;c)r(l/2 — jc) =n/ cosTU, the integrand is seen to be equal to 



r(D-i+z)r(-D/2 + i-z) /1-z 

^ — COSTTZ 



(B.113) 



up to a constant independent of z- We now perform the integral by the residue theorem, using 
the well known residues Res^=_„r(z) = ( — l)"/n!. This is seen to result precisely in the terms in 
expression (|B.111I) . Note that the pole at z = is a double pole, and this gives rise to the logarithm 
of (1 — Z)/2 and the psi-functions. □ 



C Alternative form of master integrals Mg 

For a general graph G, a corresponding Mellin-Bames type formula for the master integrals Mq 
[cf. eq. (12.301 )1 can be obtained in different ways. One method, developed in [14J, is based on 
the use of graph-polynomials and was described in the main text, leading to eq. (12.361 ). Another 
method lfT9l is to proceed by induction, integrating the vertices X^+i, . . .,Xe+v G in eq. (|2.30l) 
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one by one, and using at each step eq. (|2.32l) . We now present the result of this second method 
here, adapted somewhat from [|T9l . It yields a representation of the form 

MG(Zi,...,X£) = (47r)^(^+i)/2/^G(vv,H) n \—r^] ' (C.114) 

J'^ \<r<s<E \ ^ / 

where the w integrations are over certain contours parallel to the imaginary axis, and are absolutely 
convergent. Zy^ = Xy ■ is the Euclidean inner product on the sphere, which after analytic contin- 
uation to deSitter space becomes the Lorentzian inner product. This is of the same general form 
as the formula for the master integral given in the body of the text, (|2.36l) . but the kernel Kq^ the 
exponent Oiys, and the integration variables w, are defined differently here. To describe them, one 
first needs to introduce a somewhat elaborate notation. 

The graph G may have multiple lines going between a pair of vertices /, j. We replace such 
multiple lines by a single line, which carries a corresponding complex parameter Zij as in eq. (12.291 ). 
We call the new graph G as well for simplicity. We next choose an ordering of the integration 
vertices i = E + 1, . . . ,E -\-V which specifies the order in which the subsequent integrations of the 
corresponding X/'s are done in eq. (|2.30l) . For each z = £ + 1 ,...,£■ + V, we introduce two auxiliary 
graphs Cy, Gj, J = 1 , . . . , V as follows. The graph Cj is a 'complete graph' , meaning that each vertex 
is joined to any other vertex by one line. The vertex set of Cj is a subset of vertices from G and 
is obtained as follows. We first consider the interaction vertices {£■ + V + 1 — j, ...,£■ + V} of G, 
and call this set Vj. Then we let the vertex set 'VCj of our complete graph to be the set of all k 
such that {ik) G G for some / G Vj. The edge set "LCj consists of all pairs {ik) with i,k G I^Cj, 
because Cj is by definition a complete graph. The graph Gy is in a sense the complement of Cj. 
More precisely, the vertex set I^Gj consists of the vertices of G minus the vertices of Vj, while the 
edge set is EGj = {{ik) \ {ik) G G, i,k ^ Vj}. For each edge {ik) in Cj, we introduce a complex 
integration parameter Wjj^ G C, and we denote by w the set of all such parameters as j =\,...,V. 
Our formula for Kg is then 

Kg{z,w):^ (C.115) 
T\jT\k(^Cj^{-v^j,k{E+v-j) + 'Li(:,Cj+y Wj^i^ik) Y\jY\k(^Gj^{-Zk{E+v-j) + 'Lic,Cj+y 

Y\j^{-I.k<EGj^k{E+V-j) +Y.keCjWj,k{E+V-j) +D) 
r[ 'Lk<EGj{-Zk{E+V-j) + iE/eCy+i Wj+I4k) +I.k<ECj{Wj,k(E+V-j) -\Y.i<ECj+i Wj+\,ik) +D/2 ] 



n,;r(-z,,) 

Wj>QWi,k=U...,E: (ik)eCj ^{-Wj,ik) 



dw 



Ini 



The sums/products over j run from to V unless stated otherwise, and for j = 0, Cq is defined to 
consist of the lines {{E + V)k), where k are vertices in G connected to the vertex E + V. For those, 
we are using the notation Wo,(£+v)yt = Z[E+v)k- Our formula for a„- is 

c^rs{w):= ^J." ■ (C.116) 

;>0: Cj3(rs) 
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The integration over w is over contours parallel to the imaginary axis, such that the real part of any 
argument of a gamma-function in the numerator is > 0. Such contours do not exist for all choices 
of Zij, but e.g. for those with — £ < < for sufficiently small £. For other choices, the master 

integral is defined by analytic continuation. The vv-integrals are absolutely convergent. 

For example, for V = 1 (G the 'star graph'), the complete graph Ci consists of all edges {{ik) : 
1 < i < k < E}, and the kemel becomes 



'J) 



rizi + ...+ZE+D)r{-zi)..T{-ZE) 

(C.117) 

in accordance with eq. (12.321) . 

The representation (|C.l 141 ). and eq. (|C.l 151 ) can be proved by induction on the number V 
of interaction vertices of the graph. Each time a new interaction vertex is added, we may use 
eq. (|2.321) . The inductive process is very similar to that described in [19] . In the course of the 
inductive argument, one also proves that the integrals over w are absolutely convergent. In more 
detail, if the graph G has V vertices, it is convenient to rename the 'last' set of integration variables 
Vrs '■= wv,rsifiS = l,...,E, and to define the kemel Aq(z,v) to be the kemel Kg{ZtW), divided by 
Y[r<s r(— Vri) , and with all the integrations over j <V carried out already. The master integral 
Mg is then, after shifting the integration variables: 

MG(Xi,...,X£)=cst.G / Ag(z,v) n ^i-^rs) (^—^) (C.118) 

■^v i<r<s<E \ ^ / 

The induction hypothesis may then be formulated in terms of this kemel. They are II) that Ag(I, v) 
is analytic in v in a strip — £ < J(v„) < [where the integration contours mn in the previous 
formula], and 12), that 



< est. 



1+ E I^MI 

l<r<s<E 



e^pf-?E E^M ) ' (C.119) 



for some A^, and for sufficiently smalj^ |3(z/j)|. This condition, together with the growth estimate 

for the gamma-function, |r(-v„.)| ~ cst.|J(vr^)|"^/^+^(''")exp[-||J(v„.)|] for large p(v„)| -)-oo, 
is immediately seen to imply the absolute convergence of the integrals over v in eq. (1C.1181) for 
example if all the points X,- are in mutually spacelike position, i.e. Z„ < 1 for all r < s. More 
general configurations are discussed below. 
One has: 



Theorem 4. The master integral Mg is given by eq. (IC.l 141) . with kemel Kg(z,w) given by 
eq. (|C.l 151) . and a„ by (IC.l 161) . The multiple integrals over w are along contours parallel to 
the imaginary axis leaving the left resp. right poles of any gamma-function in the numerator to the 
left resp. right. Such contours exist e.g. for — £ < yi{zij) < 0, and the corresponding integrals over 
w are then absolutely convergent for mutually spacelike configurations of points, i.e. Z„ < 1 for all 



^Note that Mq is only needed for such Zij, as the variables Zij are related to ze by eq. (|2.29l l, and each ze is integrated 
over a contour c around 0, — 1 , —D/2 + 1 which may have arbitrarily small imaginary part. 
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r < 5, at least when \3{zij) \ < £. For general values and \'3{zij) \ < £, Mq can be analytically 

continued to a (meromorphic) function. 

Proof: The proof is by induction in V using the inductive assumptions II) and 12), and using 
at each step eq. (I2.32I ). The details are given [|T9l : a difference is that our inductive bound 12) 
is sharper than that used in [|19l . Another difference is that we give an explicit formula for Kq. 
The bound (IC.l 191 ) is seen to reproduce itself using the estimate on the gamma-function in the 
induction step, and it is satisfied for the integral (|2.321 ). □ 



In order to define Mq for more general configurations, one has to investigate at the convergence 
properties of the integrals (IC.l 141 ) or (|C.l 18t . Assume we have a complex deSitter configura- 
tion {Xi,...,Xe) in Zrs G C \ [1,°°), and let (p„ := Arg(l — Z„) G (—71, tt). In view of 
1(1 -Z„)^'" I = |1 -Z,.,|^(^'") e<P"^(^"), and of eq. (|CTT9|) and (ICTTSl) . we need to look at the 
convergence of the integral on the right side of 



|Mg(Xi,...,X£)| <cst. Yl |1 

l<r<.s<£' 
- N 

1+ E l^(vr.)l exp 



7 I sup 5H(vr. 



(C.120) 



l<r<s<E 



n ir( 

l<r<s<E 



-Vrs) \ e' 



<P/-j3(v/-i) 



The gamma functions are estimated as usual by |r(— v„) | ~ cst.| J(v„) l^^/^^^*-^"^ exp[— 1| J(v„)|] 
for large p(vri) | — )■ °°. Then, we have absolute convergence automatically if each | (p„ | <n/ 2, even 
without using the "exp" factor under this integral. This includes any real configuration (Xi , . . . , Xg) 
where all points are mutually spacelike to each other, i.e. Z„ < 1 for all r ^ s. A slightly more 
careful look at (IC.120I ) reveals that we have absolute convergence also e.g. for the following 
configuration needed in the main text. Let r < E he fixed, and let us assume that Xf is real and 
timelike w.r.t. to Xs for all s ^ r, and that the Xs are real and mutually spacelike related for 
s r. Consider instead of (Xi, . . . ,Xe) the complex configuration {X[, ...,Xr + i{r — l)ee, ...,Xe + 
i{E — l)ee), where e E MP^^ is some time-like vector and £ > small. Then for the complex 
configuration \n — (Prj| ~ cst.£, r s, and \(Pst\ ~ cst.£, r ^ sj. The bound (IC.120I) is now seen to 
give 



MGiXi,...,Xr + i{r-l)ee,...,XE + i{E-l)ee) 



<cst. Yl |1 



-M 



(C.121) 



for some M, this time using the "exp" factor to get the bound. Bounds of this kind, together with 
the analyticity of Mq, imply using standard results in distribution theory (see e.g. ch. IX of IfTSl , 
and also the proof of thm. [T]for details) that Mq is a distribution in an open neighborhood of the 
real deSitter configuration (Xi, . . . ,Xe). 
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